Local-field effects in radiatively broadened magneto-dielectric media: negative 

refraction and absorption reduction 



Jiirgen Kastel and Michael Fleischhauer 

Fachbereich Physik, Technische Universitdt Kaiserslautern, D-67663 Kaiserslautern, Germany 

Gediminas Juzcliiinas 

Institute of Theoretical Physics and Astronomy, Vilnius University, A Gostauto 12, Vilnius 01108, Lithuania 

(Dated: February 1, 2008) 

We give a microscopic derivation of the Clausius-Mossotti relations for a homogeneous and 
isotropic magneto-dielectric medium consisting of radiatively broadened atomic oscillators. To this 
end the diagram series of electromagnetic propagators is calculated exactly for an infinite bi-cubic 
lattice of dielectric and magnetic dipoles for a lattice constant small compared to the resonance 
wavelength A. Modifications of transition frequencies and linewidth of the elementary oscillators are 
taken into account in a selfconsistent way by a proper incorporation of the singular self-interaction 
terms. We show that in radiatively broadened media sufficiently close to the free-space resonance 
the real part of the index of refraction approaches the value -2 in the limit of pA 3 S> 1, where p is 
the number density of scatterers. Since at the same time the imaginary part vanishes as 1/p local 
field effects can have important consequences for realizing low-loss negative index materials. 

PACS numbers: 



INTRODUCTION 

It is well known that in dense dielectric materials the 
induced polarization P alters the field strength Ei oc act- 
ing on the constituents (i.e. the local field) compared to 
the average macroscopic field E m . Macroscopic consider- 
ations show that in systems with high symmetry such as 
a cubic lattice the two fields are related to each other ac- 
cording to E loc = E m + P/(3e ) UM- This lcads t0 thc 
well-known Clausius-Mossotti relation for the permittiv- 
ity e{uj) 



1 



pa(aj)/e 
1 - pa(oj)/(3e ) 



(1) 



where p is the density and a(uj) the polarizability of the 
oscillators. Similar arguments hold for a purely magnetic 
material except that the required densities are usu- 
ally much higher due to thc smallness of magnetic dipole 
moments and polarizabilitics. In linear response ev(u>) is 
well described by a damped-oscillator model [l[ 



a{uj) = a' + i a" = 



ato 



(2) 



The corresponding (real-valued) parameters such as the 
oscillator strength ao, thc resonance frequency and 
width, ujq and 7, are determined by the microscopic 
model. In general the linewidth 7 contains radiative 
as well as non-radiative contributions. For purely ra- 
diative interaction these parameters are strongly af- 
fected by the renormalization of energy levels and spon- 
taneous emission processes caused by the interaction 
with the vacuum electromagnetic field in the medium 
0, S E 0, IE H E3, EH- Since thc mode structure of 
the electromagnetic field inside a dense medium can be 



substantially modified compared to free space, one would 
expect that the polarizability entering cq.([T]) is different 
from that in free space. In a macroscopic approach a(uj) 
is however an input function and no conclusion can be 
drawn about possible changes due to the different struc- 
ture of the vacuum modes inside the medium. To take 
into account the modification of transition frequencies 
and radiative linewidth in a dense medium in a self- 
consistent way requires a microscopic approach. 

In the present paper we develop a microscopic ap- 
proach to local field effects in dense materials with simul- 
taneous dielectric and magnetic response using Greens- 
function techniques similar to those used by de Vries and 
Lagendijk for purely dielectric materials [l2|. To this end 
we consider an infinitely extended bi-cubic lattice of elec- 
tric and magnetic point dipoles with isotropic response 
with a lattice constant small compared to the transition 
wavelength. We however do not make use of the as- 
sumptions made in [l2l | to renormalize the singular self- 
interaction contributions to the lattice T-matrix which 
eliminated radiative contributions to linewidth and tran- 
sition frequencies altogether. We show that instead the 
self-interaction contributions can be summed to yield the 
dressed t-matrix of an isolated oscillator interacting with 
the vacuum modes of the electromagnetic field in free 
space. In this way we derive Clausius-Mossotti relations 
for general, radiatively broadened, isotropic magncto- 
dielectrica. Apart from non-radiative broadenings, the 
electric and magnetic polarizabilitics entering these equa- 
tions are shown to be exactly those of free space. We 
then show that simultaneous local-field corrections to 
electric and magnetic fields in purely radiatively broad- 
ened magneto-dielectrica have a surprising and poten- 
tially important effect: For sufficiently large densities the 
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real part of the refractive index saturates at the level of 
—2. At the same time, the imaginary part of the com- 
plex index approaches zero inversely proportional to the 
density. Thus the medium becomes transparent and left- 
handed i.e. displays a negative index of refraction with 
low absorption. 



LOCAL-FIELD EFFECTS AND 
RENORMALIZATION OF RADIATIVE 
SELF-INTERACTION IN DIELECTRIC MEDIA 

We start by developing a microscopic scattering ap- 
proach to local-field effects in dielectric media taking 
into account possible material induced modifications of 
radiative lincwidth and transition frequencies in a self- 
consistent way. To this end we consider a simple cubic 
lattice of electric point dipoles with isotropic bare polar- 
izability ab 



ab(r) 



R 



(3) 



where R denote lattice vectors. The dipoles interact with 
the quantized electromagnetic field E which obeys the 
vector Hclmholtz equation 



V x V x E(r,w) 7 E(r,w) = (i uj 2 P 



(4) 



In the weak-excitation, i.e. linear response limit, the op- 
erator of the microscopic electric polarization P has the 
form P(r) = at(r)E(r, u). Solving cq.([l]) we can deter- 
mine the (isotropic) dispersion relation k = k(u>) from 
which the permittivity e(co) can be extracted. In the 
linear response limit the solution of the quantum me- 
chanical interaction problem can most easily be obtained 
by means of Greensfunction techniques. In particular it 
is sufficient to calculate the scattering T-matrix of the 
oscillator lattice. The dispersion relation can then be 
obtained via [13|, uA, LL5| 



dctT~ 



0. 



(5) 



The scattering T-matrix obeys a linear Dyson equation 

t = v + vg (0) v + ■■■ = v + vg {0) T, (6) 

where (/^(r, r', ui) is the free-space retarded propagator 
of the electric field which is a solution to the classical 
vector Hclmholtz equation 

V x V x 0(°)(r,r',a/)- i^-0(°)( r> r', w) = 



and 



1 8(r - r') 



V(t,u) 



w 2 a b (r) 



e c z 



(7) 



(8) 



is a linear, isotropic point vertex. Note that integration 
over spatial variables was suppressed in eq.© for nota- 
tional simplicity. 

For a cubic lattice of isotropic scatterers, the series can 
be summed up to yield [l6| 



T(k,k') = - 



£ e *(k-k')R'j 1 + £ e *'R a (0) (R) j t 

(9) 

where (R) stands for Q(°> (r, r + R, loq) which due to 
the discrete translation invariance is independent on r. 
The single-particle scattering t-matrix t(to) is determined 
by the bare polarizability [12j 



(10) 



Note that g(°\0) is diagonal and isotropic. In cq.© we 
have separated the contribution of the lattice Er^o) 
from the multiple scattering events at the same oscilla- 
tor (G^ (0)). This separation is crucial since Q^°\0) is 
singular. Rather than eliminating this singularity by a 
regularization procedure as done in 12j , we note that ex- 
pression (|10p gives the single-particle scattering t-matrix 
t(w) dressed by the interaction with the vacuum field in 
free space. This quantity is experimentally observable 
and is related to the single-particle polarizability a(u>) in 
free space: 



a(u) = t(u}) -^e 



111) 



at on the other hand is not observable and thus only a 
theoretical notion. At this point other broadening mech- 
anisms can be incorporated by adding appropriate non- 
radiative decay rates -y non - rad to the polarizability a(uj) 
(fTTjl (cf. equation j2]) and discussion thereafter). 

Obviously, for the radiative part separating the sum 
J2r e lk ' R S (0) (R) into gW (0) + £ R#0 e^'V (R) does 
the trick of writing the full lattice T-matrix in terms 
of the known free space t-matrix. As a drawback we 
are left with the sum over the lattice vectors R / 0. 
Unfortunately this sum can not be evaluated exactly and 
has to be treated approximately. 

According to Poisson's summation formula 



£ f(n)= / d */( 



(12) 



fc— — OO ' 



the sum over R^0 can be expressed in terms of a real 
space integral and a sum over inverse lattice vectors K 
of the Fourier transform of the free space Greensfunction 

0P)(p) 



^e 4kR £(°)(R)=^ /drdp 

R#0 K 



S(|r|) 
(2™) 3 



,i(p+k-K)rg( ) ^ 

(13) 
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Here H(|r|) is some smooth function with 5(0) = and 
S(|r| > 0) — > 1 introduced to prevent the integral from 
touching the excluded singular point r = 0. 

In the following we restrict the discussion to lattices 
with a lattice constant much smaller than the resonant 
wavelength, i.e. ka -C 1. In this limit the lattice of 
oscillators behaves essentially as a homogeneous medium. 
Contributions from large K-vectors to the sum, which 
reflect the discreteness of the lattice, can be neglected as 
long as the singular contribution from the origin has been 
excluded. Therefore we only keep the term K = and 
assume a Gaussian cutting function E(|r| 
with S <C a. This yields 

^— ' a 3 

R#0 

1 TT 3 / 2 5 3 



1 



-r 2 /6 2 



a 3 (2tt) 3 



(14) 



where p = p/|p|. Apart from the Gaussian p-integral 
which provides a smooth cut-off in reciprocal space, 5 
can be treated as a small parameter. That allows to 
carry out the integration analytically which in leading 
order of 8 yields 



1 



y e *Rg(o) (R)RJ J_3(o) (fe) 



i 



i 



R#0 



3w 2 /c 2 



The free-space Greentensor CJ(°)(k) is given by [lj 

/ 2 \ -1 

g(o)( k ) = r^.i_|k| 3 A k 



(15) 



(16) 



with At = 1 — k ® k being a projector to directions 
orthogonal to k. 

With this we are ready to evaluate eq. ([5]) which reads 
in the limit ka <C 1 



det 



1 



1 



'-pa(uj)/e 



|kl 2 At 



0. 



(17) 

Solving eq. (|17[) for the (isotropic) dispersion k = k(u>) 
with k(u) = e(uj)Lu 2 / c 2 finally yields 

pa(u))/£ 



e{uo) = 1 + 



(18) 



1 - pa{uj)/3e 

This is the well-known Clausius-Mossotti relation where 
for purely radiatively broadened systems a(u) is the 
dressed polarizability of an isolated oscillator interacting 
with the free-space electromagnetic vacuum field. 



LOCAL-FIELD EFFECTS FOR 
MAGNETO-DIELECTRICS 

We now extend the above discussion to the case of a 
bi-cubic lattice of electric and magnetic dipole oscillators. 



The microscopic, space-dependent bare electric polariz- 
ability ab e (r) is then given by 



a be (r)=a be ^<5(r-R) = ^]T e lKr (19) 



R 



K 



and, similarly, the bare magnetic polarizability by 



a.bm{*) = OL hm 5(r- R- Ar) 



R 



a bm \ ^ e ,K(r-Ar) 

a 3 ^ 

K 



(20) 

Here R denotes again the lattice vectors and Ar the spac- 
ing between the electric and magnetic sublattices. The 
bare atomic polarizabilities ctbe and abm are assumed 
to be scalar for simplicity corresponding to an isotropic 
medium. The last expressions in eqn. (fT9|) and ([20]) give 
the bare polarizabilities in reciprocal space, with K being 
the reciprocal lattice vectors. 

Due to the simultaneous presence of electric and mag- 
netic dipole lattices we now have to solve the coupled set 
of vector Helmholtz equations for the operators of the 
electric and magnetic fields 



VxVxE T E = ioj/i Q W x M + p uj 2 P (21) 



and 



2 2 

V x V x H — ^-H = %rM. — iluV x P. (22) 



In linear response the operator of the polarization P 
and the magnetization M are proportional to the elec- 
tric and magnetic fields respectively, P(r) = af, e (r)E(r) 
and M(r) = ^ a &I „(r)H(r). 

In the following we will pursue a slightly different ap- 
proach to solve the coupled set of equations than used 
in the previous section. Taking into account the lattice 
symmetry we first write the field variables in the form 



E(r) = f dk VE(k-K 

Jl.BZ K 



)e 



i(k-K)r 



(23) 



where the dependence on frequency u> was suppressed for 
notational simplicity. The subscript denotes integration 
over the first Brillouin zone. Substituting this and the 
corresponding expression for H into (j2"l"|) - ([22"|) gives the 
Helmholtz equations in reciprocal space. After some ele- 
mentary manipulations the following closed set of equa- 
tions is derived: 
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1 

'-pa be /e 



E 



K c 



■l-|k-K| 2 A k _ 



k K 



yi(k-K') = ^y^ 



iKAi 



(k-K)x 



^ ^-1- |k-K| 2 A k _ K j 



^H(k-K')e- iK ' A ' 



-PMoQ6 m ^ ^l-|k-K| 2 A k _ K 



^H(k-K')e 



-iK'Ar 



C abe 



J K' 



E 



iKAi 



(k-K)x 



WAioOftm ^ ^1 - |k - K| 2 A k _ K j 



(24) 

^E(k-K') 



r 



(25) 



where p = 1/a 3 is the particle density. The sum in the 
brackets on the left hand sides of eqs. (|24I25[) can be 
rewritten as 



•E 



i 



k £l-|k-K| 2 A k _ K ^ 

= g(°)(0) + ^ e * kR s (o) (R), 

R#0 

where in the second line we have separated the singular 
contribution Q^(0). One recognizes that this term can 
be added to the expressions containing the bare polariz- 
abilitics in cqs. ([2"4"|) and (|23|) yielding the dressed scatter- 
ing i-matrices for isolated electric and magnetic dipoles 
interacting with the free-space vacuum field: 



LP 2 &be 

c 2 e 
u 2 Po 



S (o) (0), 



(26) 
(27) 



The sum over the Greensfunction excluding R = can be 
evaluated in a similar way as in the previous section. If 
we again assume a lattice constant a much smaller than 
the resonant wavelength, reciprocal K vectors different 
from zero can be disregarded. This leads to 



f 



PO^bm 



kx 



1 



_pt m (u>) 



UJCtbe 



C 2 Otbe 



■l-k 2 A k 
1 



E(k) = 
H(k), 



3tu 2 /c 2 
kx 



^Ottbrn 75-1 - fc 2 A k 



H(k) = 
E(k). 



(28) 



(29) 



Since we are furthermore only interested in propagating, 
i.e. transversal modes, we can further simplify the calcu- 
lation by projecting onto transversal modes using A k 



k 2 



A k E(k) 



P-QOtbm 



LOdbe 



k 2 



k x A k H(k) 



(30) 



f 



rppaa m (uj) 



k 2 



k x A k H(k) = 



C 2 Otbe 



UJflOOtbm 



-A k E(k). 



r 



(31) 



Here we have substituted the dressed single parti- 
cle ^-matrices by the free-space dressed polarizabilitics 
a e (m)M = t e{m) {uj)c 2 /u 2 en (pq 1 )- 

In order to find the dispersion k(uj) = n 2 ui 2 /c 2 we have 
to determine the solution of the secular equation of the 
linear set of eqs. (|30I31|) , which results in the condition 



1 



1 



1 



r P a e (uj)/e a £ - k 2 3i 
1 1 



(32) 



'-pp a m (LLi) 



k 2 



Solving for the refractive index of the transversal modes 
then gives n 2 = ep, where 



e = 1 



pa e (uj)/s 
1 - pa e (uj)/3so 

PPoa m (u) 
1 - pp a m (u)/3 



(33) 



(34) 



are the relative dielectric permittivity and magnetic 
permeability, respectively, both satisfying the Clausius- 
Mossotti relations. 
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Note that for longitudinal modes eqs. (|28|) and (|29|) de- 
couple. This can be seen by applying the corresponding 
projector to longitudinal waves k(g)k which leads to a dis- 
appearance of the cross-coupling terms. The dispersion 
obtained in this way gives either e = corresponding to 
electric excitons [17J, [l8( or /x = for magnetic excitons. 



NEGATIVE REFRACTION AND ABSORPTION 
REDUCTION DUE TO LOCAL FIELD EFFECTS 
IN MAGNETO-DIELECTRIC MEDIA 

It is interesting to consider the implications of the 
Clausius Mossotti relations for radiativcly broadened me- 
dia in the large density limit. Let us first consider a 
purely dielectric medium and let us assume that the po- 
larizability a e (uj) = a' e (uj) + i a"(u>) does not depend on 
the density, i.e. the medium is radiatively broadened. In 
this case one finds 



9 , . 1 Veoa'e 
-2> + 'l — ; irr 



P \a e 



(35) 



In the high-density limit and sufficiently close to reso- 
nance the response saturates at a value of —2 with an 
imaginary part that vanishes as 1/p. At this point the 
medium becomes totally opaque since the index of re- 
fraction attains an imaginary value n = iy/2 indicating 
the emergence of a stopping band. This is illustrated 
in the left column of Fig Q] for a medium composed of 
either electric or magnetic dipole oscillators. For small 
densities (p|ao|/3 = 1/3) the resonance is centered at u>o 
whereas for larger densities (p|ao|/3 = 3) the response 
shifts to smaller frequencies and is amplified. Eventu- 
ally (p|o:o|/3 = 30) the refractive index becomes almost 
purely imaginary in which case light cannot propagate 
any longer. 

This behavior changes dramatically if we consider me- 
dia with overlapping electric and magnetic resonances 
described by both an electric polarizability a e (tu) and a 
magnetic polarizability a m (uj). Independent application 
of Clausius-Mossotti local-field corrections to the permit- 
tivity and the permeability leads in the high density limit 
to 



(36) 



Thus in the spectral overlap region the real part of the 
index of refraction approaches the value —2, i.e. attains 
a constant negative value. Furthermore the imaginary 
part, responsible for absorption losses, approaches zero in 
that spectral region as 1/p. This rather peculiar behav- 
ior is illustrated in the right column of FigJTJ One clearly 
recognizes the emergence of a spectral region around the 
bare resonance frequency where the real part of the re- 
fractive index approaches —2 while the imaginary part is 
strongly suppressed. 
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FIG. 1: (color online) spectrum of the real (solid) and imagi- 
nary (dashed) part of the refractive index as well as the real 
(dotted) part of the response function(s) e and/or /iasa func- 
tion of the detuning A for a (a) pure dielectric or magnetic 
medium for p|ao|/3 at A = equal to = 1/3 (top), 3 (middle) 
and 30 (bottom) (b) magneto-dielectric medium for p|ao|/3 
at A = equal to = 1/3 (top), 3 (middle) and 30 (bottom). 



Negative refraction of light is currently on< 
most active research areas in photonics [ll| 20L 



one of the 
due 



21 



to fascinating potential applications such as superlenS' 
or electromagnetic cloaking [2 



24, m 



In 



ing [2 

recent years substantial progress has been made in re- 
alizing neg ative refraction in so-called meta-materials 
Ifil 27, 28, [2^]. These are artificial periodic structures 
of electric and magnetic dipoles with a resonance wave- 
length much larger than the lattice constant which thus 
form a quasi-homogeneous magneto-dielectric medium. 
In order to achieve a large electromagnetic response, op- 
eration close to resonance is needed which is associated 
with rather substantial losses. The elimination of these 
losses represents one of the main challenges in the field 
[30I ] . We have shown here that in a radiatively broad- 
ened medium, i.e. a medium in which density-dependent 
broadening mechanism can still be disregarded for suf- 
ficiently large densities, local field effects can provide a 
negative index of refraction and at the same time effi- 
ciently suppress absorption losses. 



SUMMARY 

In the present paper we have given a rigorous micro- 
scopic derivation of Clausius-Mossotti relations for both 
the electric and magnetic response in an isotropic, radia- 
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tively broadened magneto-dielectric medium formed by 
a simple bi-cubic lattice of electric and magnetic dipoles. 
As opposed to previous microscopic approaches we have 
taken into account possible modifications of the single- 
particle polarizabilities by the altered electromagnetic 
vacuum inside the medium in a self-consistent way. For 
a simple bi-cubic lattice it has been shown that the po- 
larizabilities entering the Clausius-Mossotti relations are 
those of single oscillators interacting with the free-space 
vacuum field. We showed that as a consequence of the 
local field corrections a radiatively broadened medium 
with overlapping electric and magnetic resonances be- 
comes lossless with a real part of the refractive index 
approaching the value —2 in the high-density limit. The 
latter could provide an interesting avenue to construct ar- 
tificial materials with negative refraction and low losses. 
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